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Mean Motions In Conditionally Periodic Separable 

Systems 1 

John P. Vinti 
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A search of the literature failed to disclose any general statement or proof of a 1 heorem 
informally current among dynamical astronomers. The present paper gives a proof of the 
theorem, which states that, in any conditionally periodic separable system the mean fre- 
quency rik of any separation coordinate g^ is equal to */&= dai/oJ &. Here a\ is the energy and 
J k is the &'th action variable. The proof is carried out for nonsingular Staeckel systems, 
so that it is applicable to any nonpolar orbit of an artificial satellite, when the potential 
leads to separability. 



1. Introduction 

Conditionally periodic separable systems are com- 
monly illustrated in works on advanced dynamics by 
the motion of a particle under the joint action of 
harmonic oscillator forces at right angles. This 
happens to be a very special case, for which each 
rectangular coordinate q k has a constant frequency, 
equal to the corresponding "fundamental frequency" 
i> k = dai/dJ k , where a x is the energy and J k the cor- 
responding action variable. 

In a more general system of this type, each gen- 
eralized coordinate q k may have a variable frequency, 
but it appears to be generally believed among dy- 
namical astronomers that the mean frequency of q k 
must be equal to v k , if the conditionally periodic 
system is separable. Needing to refer to such a 
theorem in solving a specific problem, I have 
searched the literature but have found no explicit 
statement or proof of it. The present paper is an 
attempt to furnish such a reference, with a proof 
sufficiently general to be applicable to all the sepa- 
rable problems that may arise in the gravitational 
theory of the orbit of a satellite of an oblate planet. 

Little of the analysis in this paper can claim to be 
really new. Much of the pertinent material in the 
literature, however, is discursive, relatively unavail- 
able, and expressed in notations now unfamiliar to 
most mathematical physicists. Some of it is inade- 
quate, if not incorrect, especially in the treatment of 
the periodicity of the q's as functions of the angle 
variables. Moreover, none of it seems to have been 
carried out in the shortest and most appropriate way 
to prove the theorem in question. The present 
paper attempts to give a concise and correct treat- 
ment that will serve this purpose. 

It is easy to see why such a theorem should have 
escaped formal statement and proof. Physicists 
have not been concerned with mean frequencies of 
this kind. Dynamical astronomers have been, but 
ordinarily for nonseparable systems. Until 1957 
their only separable problems were the Kepler prob- 

' This work was supported by the U.S. Air Force, through the Office of Scien- 
tific Research of the Air Research and Development Command. 



lem [l] 2 for which the coordinate frequencies are 
all equal to the constant value ^1 = ^2 = ^3, and the 
problem of two centers, which has remained a 
curiosity to be found in Charlier's famous book [2], 
but without application. Since 1957 various po- 
tentials have been suggested, by Sterne [3], by Gar- 
(inkel [4], and by the author [5], for the gravitational 
held of an oblate planet, all of which Lead to separa- 
bility and to intermediary satellite orbits with 
v\ 5^ H 5** v z . Since the solution for these orbits is 
greatly facilitated by knowledge of the mean coordi- 
nate 4 frequencies, it now becomes desirable to have 
a formal and general proof of the theorem. 

It is here convenient to discuss briefly the general 
plan of the paper, without definitions. The difficulty 
in proving the theorem arises only when the funda- 
mental frequencies v h . . ., v n are incommensurable. 
If w h . . ., w n are the angle variables, the plan is first 
to show that for a conditionally periodic system 
there exist infinitely many values of the time, with 
no upper bound, at which the orbit in itf-space passes 
arbitrarily closely to points separated from the 
initial point w x (0 ),..., w n (0) by integer intervals 
Aw k = m k , fc=l, . . ., n. This fact follows directly 
from a theorem of Dirichlet, which is easy to under- 
stand and to apply. 

To convert this result to g-space, it is necessary to 
know that q k is a single-valued, continuous, periodic 
function of the w's. To show this the author re- 
stricts considerations to nonsingular Staeckel sys- 
tems, proving that for them each q k is a single-valued 

A and each v k of 

the angle variables w h . . ., w n . With a careful use of 
the single-valuedness, the periodic property then 
follows. 

Application of these properties of the qx as func- 
tions of the w's shows that, at the values of the time 
mentioned above, the orbit in g-space then passes 
arbitrarily closely to points where each q k would have 
gone through exactly m k cycles. The proof of the 
theorem then follows. 



differentiable function of v k - 



}<IV_/1J- u k 10 

=J dq k /p k 



2 Figures in brackets indicate the literature references at the end of this paper. 
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2. Staeckel Systems 

If gi, . . ., q n and p h . . ., p n denote the generalized 
coordinates and momenta of a dynamical system of 
n degrees of freedom, the system is said to be of the 
Staeckel type [2, 6, 7, 8, 9] if the Hamiltonian H is 
given by 






q n )pl + V(q h . . ., q n ), 

At>0, (fc=l, ...,*) 



and if there exist functions (fraiqij^Mi), 
n } such that 

A k =M k ,/det(<t>ij), 
V=±,M2 k )A*, 

k = l 



^i=ii 



(i) 

• • > 

(2) 
(3) 



Mfci being the cof actor of <fr k \(q k ) hi the determinant 
det (0^). 

Conditions (2) and (3) are necessary and sufficient 
for the separability of a system with such a 
Hamiltonian. 

If we next define the domain Q of the q J s as the 
totality of real values of q h . . ., q n for which 
#1^0, fc=l, . . ., ??, we may then define a non- 
singular Staeckel system as one for which \l/ k (q k ) an d 
<t>ki(qk), i, fe=l, . . ., n, exist and are single-valued 
and for which det (</>^)^0 anywhere in Q. If we 
put 

*=(Mff<)) (4) 

for the Staeckel matrix, then <£> _1 exists and is single- 
valued anywhere in Q; in particular, 



(*-%=;!* (*=!,... ,n) 



(5) 



all exist anywhere in Q. (This restriction thus rules 
out polar orbits from consideration if the right 
ascension <p is one of the coordinates, since A z then 
becomes infinite on the polar axis.) 
The momenta p k are then given by 



pl=-2^ k (q k ) + 2 S *«(&)«<, (*=1, • 



(6) 



where the a's are separation constants, a x being the 
energy. (For satellite problems, where n=3, 2a 2 
and 2a z are usually denoted by a\ and a\.) 

The Hamilton-Jacobi function W is then given by 



n r Q k 

["[■ 



=S(±) 



■2W&)+2Efe(&)« 



]l/2 
dq k 



(7) 



where the sign is ± respectively as dq k ^0. 



3. Conditionally Periodic Staeckel Systems 

We call a Staeckel system conditionally periodic 
if each coordinate is either rotational or librational. 
A coordinate q k is rotational if: (i) it is an angle, 
(ii) with pl = F k (q k ) there exist positive real numbers 
c lk and c 2k such that c 2k ^F k (q k )^c u ^>0 for all real 
values of q k , and (Hi) if 



F k (q k +2ir) = F k (q k ). 



(8) 



Note that Ci k ^>0 rules out asymptotic motions and 
that the periodicity implied by (8) may reduce to 
simple constancy. The latter holds, e.g., when q k 
is the right ascension </> of an artificial satellite, since 
there does not exist any potential, depending on the 
right ascension </>, which both leads to separability 
and remains finite on the polar axis. 3 For such a 
rotational coordinate q k either p k y>c\ k 2 for all q k or 
p k <^ — c\' k for all q k . In either case 



*'*(&) = 



dq k lp k 



(9) 



is a single-valued function of q k , with derivative 
dv k /dq k existing and differing from zero for all values 
of q k . Thus q k is a single-valued differentiable 4 
function of v k . 

A coordinate q t is librational if there exist real 
numbers a i} b t , Cu, and C 2i and a real function 
Giiqt) such that 



with 
and 



pi=(ai-a t )(6*-«i)fl f «(2i), (io.i) 

&«£ G t (q t ) ^ C u >0 (aaqa b t ) (10.2) 
a^ 2,(0)^6,, (10.3) 



t7i(0) being the initial value of q t . If we again 
define v t by (9), then 



v i (q i )=± f" [(ii-at) (fit-tiGtiqM-wdgt, (11 



where the sign is ± accordingly as dq^Q, respec- 
tively. Then 



Vi 



J E 



G7 1/2 dE t , 



(12) 



where the uniformizing variable E t is defined by 
the equation 



2q i =a i +b i +(a i —b i ) cos E t 



(13) 



and the requirement that E t shall always increase 
as q t varies. By (12) and (10.2) v t is then a single- 
valued function of E i} with derivative existing and 



> See the tables in [9], pp. 656, 658, 660. 
4 Hereafter abbreviated to "s.v.d." 
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nonvanishing for all E u so that E t must be a s.v.d. 
functioD of Vi. By (13), however, q t is a s.v.d. func- 
tion of E i} so that, finally, q t is a s.v.d. function of 
Vi. Thus in a conditionally periodic Staeckel system 
any coordinate q k is a s.v.d. function of the corre- 
sponding v k . 

4. The VkS as Functions of the Angle-Varia- 
bles w k 

If we now let an increase of 2tt be one cycle of a 
rotational coordinate and a single round trip from 
a k to b k be one cycle of a librational coordinate, we 
may define the action and angle variables J k and 

w k by 



Jk 



and 



= (j>Pkdq k 



w k = i>W/dJ kf 



(14) 



(15) 



where W is now to be considered a function of the 
q's and the J's, rather than of the q's and the as. 
It is well known [10, 11] that J k and w k are canoni- 
cally conjugate, so that 



w k =da l /dJ k = v k . 
[f we also define the Jacobi variables Z?* bv 



(16) 



Bi^bW/bai, 

we obtain 

A = i OJ A Oat a = i 
where 

ca ki ^bJ k /dai 

■ 


(17) 

(18) 
(19) 


Increments dw u . . ., dw n then lead to 




fc=l 


(20) 


But, by (7) and (17), 




(21) 


or 

d-B«=S*M(ffiO<fe*i 


(22) 



by (21), (6), and (9). Also, by (19), (14), (6), and (9) 
=<£<l>ki(q k )dv k . (23) 



«*i 



If we now introduce the matrix <£>, the matrix 
fi=(a>fci), and the row matrices dv=(dvi, . . ., cfoj 
and dw=(dw u . . ., dw»), we find from (20) and 
(22) 

dv$=dv&. (24) 



For a nonsingular Staeckel system 5 <J>, $ l , and 12 
all exist at every point of Q, so that 



or 



dv=dwtt& 1 



(25) 
(26) 



Thus each derivative dVi/dw k exists and is single- 
valued everywhere in Q. Now the J's are all real, 
by (14). Thus, by (7) and (15), if the p's are all 
real, then Wand the w's are nil real; if some of the 
p's are nonreal, then W is nonreal and so are some 
of the w's. If all the w's are real, it ( hen follows that 
all the p's are real, else we should have a contradic- 
tion. Thus the domain Q, corresponding to the 
totality of all real values of the q's for which the 
p's are all real, also corresponds exactly to the set 
of all possible real values for all the w's. It therefore 
follows that each derivative dv t /dw k exists and is 
single-valued at any point in w-spaee. Thus each 
v k must be a s.v.d. function of w h . . ., w n . In sec- 
tion 3, however, we showed that each q k is a s.v.d. 
function of the corresponding v k . Thus for a condi- 
tionally periodic nonsingular Staeckel system each 
q k is a s.v.d. function j k {w h . . ., w n ). 

5. Periodic Properties of q k =f k (w u . . ., w n ) 
By (7) and (15) 

(27) 



&Pi 



dw k =^2 ^rrdqi. 



If now each coordinate 7, goes through an integral 
number m< of cycles, then by a familial- argument 

Aw*=S m t (b (dpi/dJ k )dqi=^jr £J m > ( VPidq_ i =m k . 

(28) 

Thus if each q k goes through exactly m k cycles, each 
w k increases by the integer m k . (Note, however, 
that such simultaneous increases are not always 
physically possible: this section is thus concerned 
only with the mathematical properties of the func- 
tions/^, . . ., w n ).) 

But we are really interested in the inverse problem 
where each w k has increased by an integer m k and 
we ask what has happened to the q's. Now the 
q's are uniquely determined by the w's, because of 
the single-valued property. In the situation of the 
preceding paragraph where each librational coor- 
dinate returns to its initial value and each rotational 
coordinate q t increases by 27rm f , each angle variable 
w k increases by m k . Since the w's determine the 
q's uniquely, this has the result that whenever Aw k = 
m k} Jfc=l, . . ., n, each librational coordinate returns 
to its initial value and each rotational coordinate 
q t increases by 2irmi. 



5 See appendix for examples. 
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Thus, in the inverse problem, whenever we are 



given Aw k ~- 



k=l, . . ., n, we find that each ro- 



tational coordinate q t must go through exactly m* 
cycles and that each librational coordinate qj must 
go through some integral number of cycles, tj, say. 
By (28), however, we then find T j =m j . Thus when- 
ever the angle variables w k are all increased by 
integer amounts Aw k =m k , each of the functions 
<Lk=fk(wi, • ■ •> w n) must go through exactly m k 
cycles. 

6. Mean Motions 

If in a time interval T the number of complete 
cycles passed through by any coordinate q k is N k , 
the corresponding mean frequency n k is, by definition 



n k = 



lim 



(W), 



(29) 



if the limit exists. We shall now prove that 
n k =v k =daildJ k , fc=l, . . ., n, for any conditionally 
periodic nonsingular Staeckel system. 

To do so, note that if p u . . ., v n are all commen- 
surable, there exist a positive v and positive in- 
tegers mi, . . ., m n such that 



Vk=mkVo, 



(k=l, 



., n), 



(30) 



where we may choose v to be the greatest common 
divisor of the v k s. Then, from (16) and (30), during 
the actual motion, 

w k =w k (0)+m k v t (k=l, . . ., n) (31) 

and in the time interval T=l/v Q we have 

Aw k =m k (k=l, . . ., n). (32) 

By section 5 each q k goes through exactly m k cycles 
in this time, so that in this case the motion is truly 
periodic, with period l/v Q . The mean frequency of 
q k is thus 

n k =m k /T=m k v =v k . (33) 



If the frequencies vi, . . . , 
mensurable, we may let 



v n are not all com- 



(34) 



and then at least one of the £'s will be irrational. 
Then by (16) and (34), during the actual motion 



f w k =w k {0)-\-^ k v l t. 



(35) 



We now use a theorem of Dirichlet [12], which states 
that if the set of real numbers £ 1? . . ., % n has at 
least one irrational member, then the system of 
inequalities 

\^-m k /P\<P~ 1 ~- (Jfc=l,...,n) (36) 
has an infinite number of integer solutions for P and the 



m's. Note that the solutions for P have no upper 
bound. 

To apply this theorem, consider only those values 
of the time interval T such that viT=P, where P 
is an integer that satisfies (36). In this time each 
w k increases from its initial value w k (0) to a final 
value given by 



w k (T)=w k (0)+P£ k , 
by (35). But by (36) 

PZ k =m k +Vk, \v k \<P n 
w k (T)=w k (0)+m k +rj k . 



so that 



(37) 



(38) 



(39) 



As viT=P takes on those larger and larger integer 
values corresponding to solutions of (36), each rj k 
approaches zero, by (38). Then, by (39), there 
exist infinitely many values of T, with no upper 
bound, at which the orbit in w-space passes arbitrarily 
closely to points where Aw k =m k , k=l, . . ., n, the 
m k s being solutions of (36). 
If the initial g's are given by 

2*(0)=/*M0), . . ., 1^(0)], (fc=l, . . , n) f (40) 

then the values of the q's at any of these times T are 
given by 



& ( T) =f k [w 1 (0) + mi + ih , 



., v)n(0) + m n +ri n ] 
(i=l, . . ., n). 



(41) 



As we let T=P/v 1 assume those larger and larger 
values already referred to, the g's then approach 
arbitrarily closelv to the values 



2*(D=/*M0)+m, 



., w n (0)+m n ] 



., n). (42) 



This conclusion follows from (38) and the single- 
valuedness and differentiability of the functions f k . 

Comparison of (40) and (42) then shows that the 
values qt(T) correspond to Aw k =m k , k=l, . . . , n, 
and are thus, by section 5, the values that would be 
reached after each q k had gone through exactly m k 
cycles. Now, by the definition (29), it follows that 
the mean frequency 



n k = 



■ lim (m k /T), 



(43) 



if the limit exists. But m k \T=v x m k \P and, as 
T^co, lim(m*/P) = ^, by (36). Thus, 



n k =i>i% k =v k , 



(44) 



by (43) and (34). 

Thus, for each coordinate q k of a conditionally 
periodic nonsingular Staeckel system, the mean 
frequency n k is equal to the corresponding funda- 
mental frequency v k =bai/bJ k . 
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7. Appendix 

For theories of satellite orbits, appropriate coor- 
dinates are spherical or oblate spheroidal. The 
corresponding Staeckel matrices and their inverses 
are, if x=r sin 6 cos <£, y = r sin 6 sin </>, z=r cos 6: 



$ = 



o 



lj) o 



Spherical 



— CSC 2 



4> 



"1 r- 2 
1 
^0 



r z esc 



cscr 



.2 tf\ 



or,if^=c[(f+l)(l-7 7 2 )] 1/2 cos</>, y=c[(?+l)(l-r) 2 )] 1/2 sin 0, s=c^: 

Oblate Spheroidal 
^^(r+1)- 1 -(? 2 +l)- 1 (r+1)" 2 " 
: C 2 r 7 2 (l-r 7 2 )- 1 (1-V)- 1 -(l-)7 2 )" 2 
1 j 

r C -I({i+l)( € I+,l)-l C -2( 1 _ |7 2)( ? 2 +|7 2)-l c -2 (? 2 +1) -l (l _ r7 2 ) -l 

-'? 2 (£ 2 +l)(£ 2 +'? 2 )- 1 ^(l-^^+r? 2 )- 1 (? 2 +) ? 2 )- 1 ^ 2 (l-r 7 2 )- 1 + r 7 2 (l+f)- 1 ] 
1 



For spherical coordinates <]> is most easily written 
down from the expressions for p 2 , p 2 ,, and p^ in 1 he 
Kepler problem [1], with replacement of a\ and a| 
by 2a 2 and 2a 3 . For oblate spheroidal coordinates 
<£ may be found by comparing eqs (53) and (59. 1 ) of 
[5] with eq (6) of the present paper. 

Note that 3> or $ -1 could fail to exist only when 
sin 0=0 or when 77^= 1 . This could happen only 
when the satellite goes over a pole and thus only in 
a polar orbit. Such a singularity in a polar orbit, 
however, is to be expected, since nonsingularity of 
a Staeckel system leads to the g's being differentia Me 
functions of the w's and thus of the time. In a 
polar orbit, on the other hand, the right ascension 
(/3=e</> is a discontinuous function of time, being 
constant except at polar crossings, where it changes 
by 7T. 



(Paper 65B2-52) 
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